Measurement of the Michel Parameter p in Muon Decay 
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The TWIST Collaboration has measured the Michel parameter p in normal muon decay, /i^ — > 
e'^Vei'fj,. In the Standard Model, p = 3/4. Deviations from this value require mixing of left- and 
right-handed muon and electron couplings in the muon-decay Lagrangian. We find p — 0.75080 ± 
0.00044(stat.) ± 0.00093(syst.) ± 0.00023, where the last uncertainty represents the dependence of 
p on the Michel parameter rj. This result sets new limits on the Wl — Wr mixing angle in left-right 
symmetric models. 

PACS numbers: 13.35.Bv, 14.60.Ef, 12.60.Cn 



Normal muon decay, /i+ — > e^Vei^p,, is an excellent 
laboratory to test the space-time structure of the weak 
interaction. The energy and angular distributions of the 
positrons emitted in the decay of polarized muons can 
be described in terms of four parameters - p, r/, ^, and 
5 - commonly referred to as the Michel parameters. Ne- 
glecting the electron and neutrino masses and radiative 
corrections, the differential decay rate for_positive muon 
decay is given in terms of p, ^, and S by 
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where is the polarization of the muon, x is the outgo- 
ing positron energy as a fraction of the maximum possible 
value, and 9 is the angle between the muon polarization 
axis and the positron decay direction. The fourth decay 
parameter, 7], contributes to the angle-independent part 
of the distribution if one includes the finite electron mass. 
In this Letter, the TWIST Collaboration reports a new 
measurement of the Michel parameter p. A concurrent 
measurement of the parameter 6 is described in Ref. [3 . 

The current accepted value of p, 0.7518 ± 0.0026 
is consistent with the Standard Model expectation, p = 
3/4. Any deviation from 3/4 would imply the muon- 
decay Lagrangian includes scalar, vector, or tensor cou- 




plings between left-handed muons and right-handed elec- 
trons or vice versa 0]. For example, in left-right symmet- 
ric models, the Wl — mixing angle ( is given by 
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Unlike many other limits on right-handed currents, this 
is independent of the form of the right-handed CKM ma- 
trix. Recently, p has also been related to loop corrections 
to the neutrino mass matrix Jj. For a review of muon 
decay within the Standard Model, see Ref. 

TWIST utihzes the M13 beam Hue at TRIUMF to 
transport beams of 29.6 MeV/c surface muons from pion 
decay-at-rest (P^ ^ —1) or 32.8 MeV/c cloud muons from 
pion decay-in-flight (P^ +0.25) into the TWIST spec- 
trometer. The TWIST spectrometer consists of an ar- 
ray of very thin, high precision planar wire chambers lo- 
cated within a 2-T magnetic field oriented along the beam 
direction. The spectrometer includes 44 drift chamber 
(DC) planes operated with DME gas and 12 multi-wire 
proportional chamber (PC) planes operated with CF4- 
isobutane (80-20) gas, all within a He-N2 (97-3) enclo- 
sure. The wire planes are symmetrically located up- 
stream and downstream of a 125-/im thick Mylar stop- 
ping target, with lO/g" pm of graphite painted on each 
surface. A detailed description of the TWIST spectrom- 
eter is given in Ref. [^. 



2 



After muons enter the magnetic field, tliey pass 
through a thin plastic scintillator that provides the event 
trigger. They then pass through the detector planes un- 
til they stop in the target. Decay positrons follow helical 
trajectories through the DCs and PCs, permitting their 
momenta and decay angles to be measured precisely. For 
each event, all DC and PC hits within an interval from 6 
fis before until 10 /is after the trigger time are recorded. 

During off-line analysis, the PC and DC hits are ex- 
amined to identify events in which the muon stopped in 
the target, then decayed at least 1.05 fis, and no more 
than 9 /zs, later. The delay insures that the PC and DC 
hits associated with the muon and decay positron do not 
overlap. Events are rejected if a second muon enters the 
spectrometer, or if a beam positron passes through the 
spectrometer within 1.05 /iS of either the muon arrival or 
decay time. Additional cuts include the muon flight time 
through the M13 beam line and a requirement that the 
muon stopping location be within 2.5 cm of the detec- 
tor axis. All events that pass these cuts are analyzed to 
reconstruct the decay positron kinematics. 

After track fitting, 2.3% of the events contain addi- 
tional tracks in coincidence with the decay. Extra tracks 
can arise from beam particles that are not resolved in 
time, events that scatter within the detector leading to 
two reconstructed track segments, and events that in- 
clude delta rays or decay positrons that backscatter from 
material outside the detector volume. Two algorithms 
have been developed to select among the choices in multi- 
track events. They also impose different constraints on 
events that scatter within the detector when only one 
track segment is reconstructed. All events have been an- 
alyzed using both algorithms. 

To extract the Michel parameters, the measured 
positron spectrum is compared to that predicted by a 
detailed Monte Carlo (MC) simulation. The MC uses 
GEANT ^3 to simulate particle interactions and a 
model based on GARFIELD to simulate wire cham- 
ber responses. The MC decay generator includes the ef- 
fects of electron mass, plus first-order and many higher- 
order radiative corrections 0. It also includes beam 
positrons and additional muons in the simulated events 
according to their observed rates in the data. The out- 
put from the MC is digitized and processed by the same 
analysis codes that are used for real events. 

The data reported here include a total of 6 x 10^ muon 
decay events that were recorded during Fall, 2002. Six- 
teen independent data sets were taken to explore the sen- 
sitivity of the spectrometer and analysis to a broad range 
of systematic effects. A typical data set included 3 x 10^ 
events, sufficient to determine p with a statistical preci- 
sion of ~ 0.0007. In addition, special runs were taken to 
provide data to validate aspects of the simulation that are 
difficult to test with the muon decay spectrum. Five data 
sets were taken under conditions that permit a reliable 
determination of p. Four sets were taken with a surface 
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FIG. 1: (color online) Measured positron momentum spectra 
(solid lines) near the end point are compared to Monte Carlo 
simulations (dotted lines). The curves show surface muon 
set B for -0.74 < cos 61 < -0.70, cloud muons for -1-0.70 < 
cos 6* < -1-0.74, cloud muons for -0.74 < cos6l < -0.70, and 
surface muon set B for -1-0.70 < cos9 < -1-0.74. 



muon beam. Sets A and B were obtained six weeks apart 
at a magnetic field of 2.00 T; the other sets were taken at 
1.96 T and 2.04 T. The fifth data set was taken at 2.00 
T with a cloud muon beam. Depolarizing interactions in 
the target reduced the average muon polarization at the 
time of decay to ^ 90% of the incident polarization. 

The only discrete feature in the muon decay spectrum 
is the end point. Figure ^ shows comparisons of the 
measured spectra near the end point to MC simulations. 
The observed end point falls below the kinematic limit 
of 52.828 MeV/c due to positron energy loss in the tar- 
get and detector materials. Fits to spectra in the region 
p > 52 MeV/c for various decay angles show that the 
energy loss follows the form, AE{9) = —a/\ cos9\, with 
a a constant as expected for the planar geometry of the 
TWIST detectors ^ takes on different values for 

upstream and downstream decays when the muon stop- 
ping distribution is not centered in the target. We use 
E = E„ieas + a/I cos 0| to correct both data and MC 
events for the average positron energy loss. 

Effects that distort the reconstructed positron mo- 
menta will lead to systematic errors in the extracted value 
of p if they are not simulated accurately by the MC. To 
test the simulation of energy loss for positron momenta 
well below 52.83 MeV/c, events were recorded in which 
a muon came to rest at the far upstream end of the de- 
tector. Positrons from muon decays in the downstream 
direction first spiral through the upstream half of the de- 
tector, then pass through the target, and finally spiral 
through the downstream half of the detector. Figure El 
shows the difference between the reconstructed positron 
momenta in the two halves, which measures the energy 
loss in the target and detector materials. The MC agrees 
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FIG. 2: (color online) Average energy loss of decay positrons 
through the target and detector materials as a function of 
momentum for data (closed circles) and Monte Carlo (open 
circles), measured as described in the text. 



very well with the data. Similar comparisons verify the 
MC simulation of positron multiple scattering and hard 
interaction rates. 

U p = Ph + A/9 and rj = rjH + A77, then the angle- 
integrated muon decay spectrum can be written as: 

Nix) = Ns{x,ph,7]h) + ApNApix) + A7/7Va„(.t). (3) 

The spectrum is linear in the decay parameters, so this 
expansion is exact. It can also be generalized to include 
the angular dependence. This is the basis for the blind 
analysis. The measured momentum-angle spectrum is 
fitted to the sum of a MC 'standard' spectrum Ns pro- 
duced with unknown Michel parameters pn, Vh, ^h, Sh 
and additional 'derivative' MC distributions Nap, Na^, 
and Na^s, with Ap, A^, and A^i5 as the fitting param- 
eters. The hidden Michel parameters associated with 
Ns are revealed only after all data analysis has been 
completed. Finally, a refit is performed to shift rj from 
the hidden value r]H to the accepted value. The fidu- 
cial region adopted for this analysis requires p < 50 
MeV/c, Ip^I > 13.7 MeV/c, pr < 38.5 MeV/c, and 
0.50 < I cos 6*1 < 0.84. 

Figure|3Ja) shows the momentum spectrum from set B 
in the angular range 0.70 < | cos 6*1 < 0.84. The probabil- 
ity for reconstructing muon decays is very high, as shown 
in Figs. Olb) and (d). Thus, higher momentum decays 
that undergo hard interactions and are reconstructed at 
lower momenta can lead to an apparent reconstruction 
probability above unity. Figures l^c) and (e) show the 
residuals of the fit of the decay spectrum from set B for 
two different angular regions. Similar fits have been per- 
formed to the other data sets, yielding the results shown 
in Table n 

The 11 additional data sets have been combined with 
further MC simulations to estimate the systematic un- 
certainties shown in Tables HI and UTI The largest effects 
arise from time- variations of the cathode foil locations 9] 
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FIG. 3: (color online) Panel (a) shows the muon decay spec- 
trum (solid curve) from surface muon set B as a function of 
momentum, for events with 0.70 < |cos^| < 0.84, as well as 
the events within this angular region that pass the fiducial 
constraints (dot-dashed curve). Panels (b) and (d) show the 
probability for reconstructing muon decays for two angular 
ranges, as calculated by the Monte Carlo. Panels (c) and (e) 
show the residuals for the same angular ranges from the fit 
of set B to the Monte Carlo 'standard' spectrum plus deriva- 
tives. 



and the density of the DME gas, which change the drift 
velocities and influence the DC efficiencies far from the 
sense wires. These parameters were monitored through- 
out the data taking, but only average values were used 
in the analysis. Special data sets and MC simulations 
that amplify these effects have been used to estimate 
their uncertainties for p. Other important effects arise 
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TABLE I: Results from the fits to the various data sets us- 
ing one of the track selection algorithms. Each fit has 1887 
degrees of freedom. Only statistical and set-dependent sys- 
tematic uncertainties are shown. 



Data Set p 




Set A 0.75134 ± 0.00083 ± 0.00053 


1814 


Set B 0.74937 ± 0.00066 ± 0.00053 


1965 


1.96 T 0.75027 ± 0.00065 ± 0.00055 


1951 


2.04 T 0.75248 ± 0.00070 ± 0.00060 


1804 


Cloud 0.75157 ± 0.00076 ± 0.00053 


1993 


TABLE II: Contributions to the systematic uncertainty in 
p. Average values are given for those denoted (ave), which 
are considered set dependent when performing the weighted 
average of the data sets. 


Effect 


Uncertainty 


Chamber response (ave) 


±0.00051 


Stopping target thickness 


±0.00049 


Positron interactions 


± 0.00046 


Spectrometer alignment 


±0.00022 


Momentum calibration (ave) 


±0.00020 


Theoretical radiative corrections 


± 0.00020 



Track selection algorithm ±0.00011 
Muon beam stability (ave) ± 0.00004 



from the uncertainty in the thickness of the graphite lay- 
ers on the Mylar target 9] and from uncertainties in the 
GEANT treatment of positron interactions that lead to 
the spectrum distortions seen in Fig. |31 Upper limits on 
these uncertainties were obtained from the data where 
the muons stopped far upstream. Smaller effects arise 
from the uncertainty in the orientation of the detector 
with respect to the magnetic field, the statistical uncer- 
tainties in the end point fits that are used in the mo- 
mentum calibrations, the theoretical uncertainty in the 
higher-order radiative corrections, and time- variations in 
the average muon stopping location within the target. 

We treat the chamber response, momentum calibra- 
tion, and beam stability uncertainties as set dependent 
when computing the average of the data sets since the 
underlying causes fluctuated in time. This gives p ~ 
0.75091 ± 0.00032(stat.), with = 7.5 for 4 degrees of 
freedom. The alternative track selection technique gives 
p = 0.75069. We average these results as our best es- 
timate of p, and include half the difference in the sys- 
tematic uncertainty. We then rescale the statistical un- 
certainty by \/xi and, to be conservative, consider the 
systematic uncertainties as common to the five data sets. 

We find p = 0.75080 ± 0.00044(stat.) ± 0.00093(syst.) 
± 0.00023, consistent with the Standard Model expecta- 
tion p = 3/4. This result assumes that rj is given by the 



accepted value 0|, i] = — 0.007 ± 0.013; the third uncer- 
tainty represents the change in p when 77 changes within 
its uncertainty. Within left-right symmetric models, this 
result sets a new upper limit, |^| < 0.030 (90% c.L), on 
the Wl — Wr mixing angle. 
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